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1 Introduction 

Some man-made landscapes, such as plots of land may be viewed as T- 
tessellations, that is a tessellations of a subset of the plane, where all vertices 
are degree three and with one flat angle. Adamczyk and Kieu have developed 
a random model aiming at representing such landscapes. This article proves 
that their model, and extensions thereof, are well-defined. 

I shall first shamelessly summarise their arguments for defining this new 
model. Models with T-vertices had already been studied [Mackisack and 
Miles, 1996], but they are usually built starting from point seeds from which 
segments grow until they are blocked: this situation does not really corre- 
spond to cultivated fields, and if the segments all grow at the same speed, 
the parcels cannot have different scales. 

At another extreme, Arak et al. [1993] have introduced a very general 
model for random graphs, directly defined by their measure on the set of 
graphs. The measure depends on an energy function that can be specified 
to yield T-tessellations, as detailed by Mackisack and Miles [2002]. These 
graphs have very nice mathematical properties. As a result, they may be 
sampled exactly, without resorting to Metropolis algorithm. However these 
nice properties come at a price: the tessellations are necessarily very random. 
For example, the intersection with any line is a Poisson process. It is thus 
hard to generate long-range spatial correlations. 

Adamczyk and Kieu define their model as a non-normalised density with 
respect to the model of Arak et al. [1993]. Different energies allow different 
kinds of landscapes. We may for example require that all parcels have a 
similar area, or penalise sharp angles. The price to pay is harder sampling, 
requiring Monte-Carlo Markov chain algorithms. 
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From a theoretical point of view, Adamczyk and Kieu have not managed 
to prove that their measure was finite, which is necessary for the model to be 
well-defined. This article focuses on proving that we have a true probability 
measure. 

In section 2, we describe the mathematical structure of the models, make 
some general remarks on the structure of T-tessellations, introduce the no- 
tations and our strategy. Mainly, it is enough to bound precisely enough 
the weight of all T-tessellations on a given configuration of lines. This can 
be achieved by bounding the number of different T-tessellations on those 
lines. We then count. To do that, we devise algorithms that can rebuild 
the tessellation from their input, and count the number of different inputs 
they can have. The algorithms and the final theorem are detailed in section 
3. Formal proofs of correctness of the algorithms are all delayed to the ap- 
pendix A. Finally we discuss the limitations of our proof and hint at possible 
ameliorations in section 4. 



2 Notations and strategy 

Let us first recall some definitions related to Poisson line processes on the 
plane R^. 

We may define a line D by an angle a e [0, 7r[ and a distance p G M. Let 
us consider the origin O = (0, 0) and the point P of radial coordinates (p, a). 
Then the line D{p,a) is the line orthogonal at point P to the line (OP). 
This parametrisation is usual since the lines corresponding to a homogeneous 
Poisson process on [0,7r[xR yield a stationary and isotropic process. 

We now consider a domain W in the plane. To avoid many complications 
thereafter, we require it to be convex. We denote by Xw the measure on the 
Poisson line process restricted to W, that is the image of the Poisson point 
process on D^/ C [0,7r[xR where Dyy corresponds to the set of lines in the 
plane that intersect W. The measure of this process is denoted A]^ where the 
intensity r is the expected number of lines in the process. It may be seen as 
a measure on the set Cw of sets of line on W. In particular, if we denote by 
LJy the random set of lines, its cardinal is a Poisson variable with parameter 
r: ^ 

P[#L^^ = A:]=exp(-r)^. 

A finite polygonal tessellation of is a finite partition of W into cells 
such that the border of each cell is an union of a finite number of segments 
and a part of the border of W itself. A vertex of a tessellation is T-shaped if it 
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is of degree three and two of the incident edges are ahgned. A T-tessellation 
is a polygonal tessellation such that : 

• All its inner vertices are T-shaped. 

• No two segments are aligned. 

We denote by Tw the set of T-tessellations on W. We denote by L(T), 
or simply L the set of lines that support the segments of T. 

The measures defined by Adamczyk and Kieu, which include the partic- 
ularization of the model by Arak et al. [1993] to T-tessellations, are of the 
form: 

d-fw{T)=exp{-H{T))dX^^{L{T)), 

where the energy H(T) > Cjj^L{T) is bounded from below by a linear func- 
tion of the number of lines in the tessellation ^L{T). As a remark, scaling 
W is equivalent to scaling r^. 

We aim at proving that the corresponding partition function is finite. We 
may bound it by: 

exp{-H{T))dy^{L{T)) 

Tw 



< 



Tn 



sup exp{—H(T')) 

T':L{T')=L(T) 



< I exp(-C#L(T))#{T:L(T) = L}dA^^(L) 

< J^exp(-Cfc)Ar(fc)exp(-r) — , (1) 

fc=0 

where we have used the Poisson distribution on the number of lines in and 
denoted by M{k) the maximum number of ways to obtain a T-tessellation as 
segments on k given lines, whatever the position of those lines: 

M{k) = sup # {T : L{T) = L} 

L:#L=k 

Our strategy consists in counting these T-tessellations on fixed lines, or 
rather bound their number from above. 

Since k\ > (k/e)'', the sum in the partition function (1) converges if, for 
any real a: 

Afik) = o(Pa^). (2) 
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This implies that k-th powers a'^ are free when we are counting the tessel- 
lations, as long as there is a bounded number of such factors. In particular, 
we may label each line with "red" or "blue" {2^ possibilities), or more gener- 
ally with any properties among a finite set. We may also split k indiscernible 
objects among the k lines (less than 4^ possibilities). 

To bound this number N'ik), we shall use descriptions of the tessellation, 
essentially labels on lines. These descriptions, or labelling schemes will be 
shown to characterise the tessellation, given the lines L, by rebuilding T 
with an algorithm using that description as input. A bound on the number 
of different descriptions then yields a bound on ^/{k). Many of the further 
notations are also thought to be easy of use within an algorithm. 

Let us have a closer look at T-tessellations. 

A T-tessellation T is built on a set of lines L that support its segments. We 
shall write s{l) for the segment supported by the line I G L, and conversely 
l{s) for the line that supports s. 

The endpoints of those segments can only be an intersection with another 
line, or with the border of W. So that, knowing L, the only places where 
something can happen are those intersections of lines and possibly border. 
We call those intersections events. A generic event is noted e. When specify- 
ing the event, through the lines that intersect, we write e(/,m), for /,m G L. 
Naturally the event e(/,m) is the same as c(m, /). A special case is when a 
line intersects the border B of W. Since W is convex, this happens exactly 
twice, so that t{B, I) and e(/, B) are different. Conventionally c{B, I) < e(/, B) 
for the order we define now. We write S for the set of events. 

Given k lines, we can find their k{k — 1)^ crossing points, as well as their 
2k crossing points with the border. They are almost surely distinct. We 
choose an axis along which each crossing has a different coordinate. This 
axis will be called the time axis, or indifferently abscissas axis. The cor- 
responding coordinates are called either times or abscissas. We order the 
events according to time. The ordered list of events will be denoted S°. The 
reverse-ordered list will be denoted S^. We shall use the usual vocabulary 
associated to time, such as saying that a point (or an event defined at that 
point) happens before another if its abscissa is smaller. We also use left and 
right for smaller and larger times. 

Since the crossings all happen at different times, the extremal points of 
each segment of the tessellation have different times. We say that the segment 
is born at its extremity with lower time, and dies at the other. We denote 
by (/) and xj(/) the times of birth and death of the segment s{l). The line 
crossed at the birth is the parent of the segment's line, and the segment's 
line is its child. The line crossed at death is the killer of the segment's line. 
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The segment's line is its victim. 

These relations thus give us two trees, the tree of births and the tree 
of deaths. Both have k + 1 nodes, labelled by the border and the k lines, 
and both are rooted at the border. For the tree of deaths, parent and child 
correspond to killer and victim. 

The tree of births encode all information about births, that is on leftmost 
extremities of segment xf{l). Symmetrically, the tree of kills encode all 
information about deaths, that is on rightmost extremities of segments a;J(/). 
So that rebuilding the two trees is equivalent to rebuilding the tessellation. 

Since the segments describe the tessellation, and for algorithmic purposes, 
we now think of a T-tessellation on a set of lines L as a couple of functions 
T = [xj , xj) with : L U i3 — )■ M. The border is added in the domain for 
the algorithms. 

Conventionally, to ease the writing of the algorithms, we now require that: 

• W is contained in the band of abscissas (0, 1). 

• The border is "always alive": xJlB) = and x'^{B) = 1. 

Even when they follow these requirements, all such couples of functions 
are not a tessellation, let alone a T-tessellation. We shall dub prototessella- 
tion any such couple. The notion will be mainly useful for initialisation of 
the algorithms. 

A T-tessellation is a prototessellation P with the following three proper- 
ties: 

• Segments do not cross: 

V/, m G L : -1 1 [xf (/), xf (m) < e(/, m) < x^(/), xf (m)] 

or [xf (/) = xf (m)] or [xf (/) = x^(m)] or [x^{l) = x^(m)] } . 

(3) 

• Segments are born on the inside of another segment, or on the border: 

If xf (/) = e(/,m), then xf (m) < e(/,m) < x^(m). (4) 

• Segments die on the inside of another segment, or on the border: 

If x^(/) = e(/,m), then xf (m) < e(/,m) < xf (m). (5) 

A prototessellation where segments do not cross (3) is a pretessellation. 
We deal with such objects within the algorithm, in some cases as output. 
We shall usually write P for all prototessellation, and xf and xf for the 
corresponding times of birth and death. 



5 



3 Algorithms and Result 



3.1 Preliminary algorithm 

If we know the tree of births alone, we can almost rebuild the tessellation. 
We only need the number of murders of each line, which is essentially free 
(4^^). Hence counting the number of tessellations in the worst case N'{k) 
is essentially equivalent to counting the highest possible number of trees of 
births on k lines. 

We start with a labelling scheme that only gives existence of Z for low 
intensities (A < (4e)~^ with non- negative energy H). However it is a basis 
of our final labelling scheme, and the proof of its efficiency introduces ideas 
that we shall use again, while staying in an easier context. 

A first way to label each line so that we may rebuild the tessellation is 
the following: 

• Give for each line its parent. That yields possibilities. 

• Give for each line its number of murders. That yields {^^") < possi- 
bilities. 

The number of murders of a line / is defined as 

M^il) = # {m : xj(m) = e(m, /) } . (6) 

Conventionally, we may set M^{B) to infinity. Though it is not necessary, it 
avoids to keep track of it in the algorithm. 

We may now rebuild the tessellation with Algorithm 1. 

Informally, we move along the abscissas axis, while prolongating the seg- 
ments that are alive. We know when each segment is born, so we add them 
then. When two segments cross, we look at their remaining number of mur- 
ders. One of the two must be zero. The corresponding segment is killed. The 
other segment's number of murders is decreased by one. When a segment 
hits the border, it is also killed. 

Lemma 3.1. Algorithm 1 yields T . 

Proof in appendix. 

The lemma yields Mik) < (ik)^. Putting that back into bound (1) would 
yield a convergent series if A < 1/ (4e) and H > 0. Compare with calculations 
in section 3.3 for details. 
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Algorithm 1 Rebuild from tree of birtlis and number of murders. 



Input: Tlie set L of lines of the tessellation, a prototessellation [xbiXd] 
such that Xft = x^, a murder function M : L U S -> N such that M = 
the number of murders in the real tessellation (6), the ordered set S° of events. 



1: for all / G L do 

2: Xd{l) ^1 > Death time temporarily set to maximum 

3: end for > End of initialisation 

4: for e G do > Consider potential events timewise 

5: /i,/2^/(e) 

6: if Xb{li) , Xb{l2) < e < Xd{li),Xd{l2) then > Do the lines cross? 

7: if M(/i) = then > Which line is killed? 

8: Xd{h) ^ e > Kill it 

9: M^k) ^ Mik) - 1 > Count that k killed it 

10: else 

11: Xdih) ^ e 

12: M(/i) ^ M{k) - 1 

13: end if 

14: end if 

15: end for 

return P 
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3.2 Main algorithm 



The previous labelling scheme still uses too much info for proving existence 
of our measures. Namely specifying the whole tree of births dooms the effort. 

Next algorithm rebuilds T while knowing only part of the times of birth. 
We call orphan a line whose parent is not specified in the labelling scheme. 
The price to pay is higher complexity: instead of one pass on events, we have 
to loop back and forth in time, prolongating orphan segments to their birth, 
and cutting too old segments until stabilisation. 

Let us be precise. The algorithm is supposed to take the following input 
data on the lines. This list will be referred to as requirements. 

1. Label each line as "leaf" or "not leaf". Leaves will not be allowed to 
have any child. 

2. For each line, its parent is given, and/or it is specified as the parent of 
its first child. 

3. Each line is given a number of virtual murders, whose properties we 
make precise after formulating the algorithm. 

4. Each line is given a number of other children, also made precise below. 

The first set of labels yields a 2'' factor. The third and fourth sets of 
labels each yield as usual less than 4*^ possibilities, as we shall see. We plan 
on showing that we can rebuild the tessellation while specifying the parent 
of few enough lines to satisfy bound (2). 

We shall denote by Uq the set of (orphan) lines without an explicit parent. 

Informally, we first initialise the algorithm by moving along the time 
axis, while prolongating the segments that are deemed alive. Namely we 
prolongate a segment when it is born, if we know the time, and prolongate 
it when it has a known child otherwise. When two segments cross, we stop 
each one if its virtual number of murders is zero. Both may be stopped at 
the same time, and at least one must be. If a segment is not stopped, its 
virtual number of murders is decreased by one. When a segment hits the 
border, it is killed. End of initialisation. 

Now we loop. Each loop consists of one pass backwards in time, and a 
cutting phase. 

The pass backwards in time is given in Algorithm 3. During the pass back- 
wards in time, we extend the segments whose parent we do not know. Since 
these segments are extended before their first child (requirement 2) in the 
real tessellation T, they never cross each other in the process. When such a 
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backwards segment hits another segment, we decrease that segment's number 
of other children. 

The cutting phase is given in Algorithm 4. During the cutting phase, 
we cut the segments whose number of other children n is negative. More 
precisely, their rightmost point is now their (— n)-th rightmost crossing with 
an orphan segment. Their number of other children is reset to zero. The 
consequence of this operation is that (— n) orphan segments will be further 
extended in the next loop. 

End of loop. Stop when all numbers of other children are zero. End of 
algorithm. 

Let us give some intuition. The algorithm hinges on a few facts. 
After initialisation, P is a pretessellation. 

Most importantly, the following property is true at the end of each loop: 

Late Events Property. Deaths occur after the time they occur in the real 
tessellation. Same for births. 

In formula, if we denote xf and the times of birth and death at the 
end of a loop, then for all lines I ^ L: 

x^l) > xUl), (7) 
x^Al) > 4(0- (8) 

Finally, within a loop, birth and death times are always non-increasing. 
That is, if we write xl and xj^ for times of birth and death after i-th change, 
then 

xl{l)>xl-'\l), 
x\{l)>xf\l) 

for all lines I E L. 

We may now state how the number of virtual murders and other children 
will be defined as input. A first consequence of these choices will be that the 
algorithm actually runs. 

For other children, we merely use the number in the real tessellation, that 
is the number of orphan segments born from the segment. In formula: 

0^(/) = # {m G ?7o : xlim) = c{l, m)} . (9) 

The definition also holds if / is the border B. Here again, we have the crude 
bound of k other children, and ways to split them. 

For virtual murders, we define them by running the initialisation phase 
without them. When two segments cross, we check on the real tessellation 
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Algorithm 2 Rebuild from final labelling scheme 



Input: The lines L, the ordered and reverse-ordered list of events E° 
and i^*", a subset U^) ^ L oi orphan lines satisfying requirement 2, a 
prototessellation P = {xb,Xd) such that Xb{l) = xj{l) for all non-orphan line 
/ G L\Uo, a function "virtual murders" V : LU B ^ N defined in (13) , and 
a function : LUB eN giving the number of orphan children a line has (9). 



1: U ^Uo 

2: for all / G L do 

3: Xd{l) ^1 > Death time set to maximum for now 

4: end for 

5: for all / G [/ do 

6: Xf,{l) ^1 > Birth of orphans set to maximum for now 

7: end for > End of "preinitialisation" 

8: for all e G ^° do 

9: /i,/2^/(e) 

10: if [xbih) = e] A [xbih) > e] then > Is /i the first child of ^2? 

11: Xbih) > Temporary maximum birth time for I2 

12: else if [xbih) = e] A [xbih] > e] then 
13: Xb{h) ^ e 

14: else if Xb{li) , Xbih) < 2 < Xd{li),Xd{l2) then > Do the lines cross? 
15: if V{h) = then > Is h "virtual-killed"? 

16: Xd{h) ^2 » Death time to new maximum 

17: if Vik) = then > Is k "virtual-killed"? 

18: Xdih) t> Death time to new maximum 

19: else 

20: Vik) ^ Vik) - 1 > Count that k "virtual-killed" h 

21: end if 

22: else > In that case I2 is "virtual-killed" 

23: V{li) ^ V{li) - 1 > Count that h "virtual-killed" k 

24: a;(i(^2) ^ 2 > Death time to new maximum 

25: end if 

26: end if 

27: end for > End of initialisation 

28: repeat 

29: P ^ Parent _seek{L, P, S°, U) > Extend backwards 

30: Cuts, U, P -ir- Cutting{L, P, E^, Uq, O^) > Cut too long segments 

31: until Cuts = 
return P 
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Algorithm 3 Parent-seeking loop 



Input: Tlie lines L, a pretessellation P = {xb,Xd), the ordered sequence of 
events S°, a subset t/ G L of lines whose parent is not currently known. 

1: for all e G do > Reverse timewise 

2: h,k^Kc) 

3: if [h eU] A [xbih] > e] A [xbik) < e < Xdik)] then 
4: Xb{li) ^ t > Extend li backwards 

5: U ^ U — li > /i seen as born on I2, for now 

6: end if 

7: if [I2 eU]A [xbih) > e] A [xb{li) < c < Xd{h)] then 

8: Xbih) ^2 > Same as above, li and ^2 switched 

9: U ^U-k 

10: end if 

11: end for 

return P 



whether a segment needs to be prolongated, and increase its virtual murders 
counter accordingly. We use the number of virtual murders we get in the 
end. An immediate consequence is that there are at most k virtual murders, 
so that there are less than 4^ ways to split them. 

Rather than copying the initialisation phase of Algorithm 2 and replac- 
ing the statements on V, let's give another characterisation of these virtual 
murders. The pretessellation P{Uq) after initialisation has the following prop- 
erties (we shall not prove it but merely use it for definition): 

xf(^o)(/)=a;^(/) if/^t/o- 

(10) 

xl^^°\l) =inf {e(/,m) : e(/,m) = x^m)} if / G Uq. 

(11) 

xJ(^°)(Z) = inf |e(/,m) : xj^^°^(m) > e(/,m) > xj(/), xf ^^°^(m)} for aU / G L. 

(12) 

Notice that the times of death are well-defined. Then the number of virtual 
murders is simply the number of kills in P[Uq), given that simultaneous 
deaths do not count: 

V{1) = # |m G L : e(/,m) = x^^^'\m) < xj^^"^(/)} . (13) 
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Algorithm 4 Cutting loop 

Input: Tlie lines L, a pretessellation P = {xb,Xd), the reverse-ordered 
sequence of events S^', a subset Uq E L oi orphan lines, a function 
: L U B ^ N giving the number of orphan children a line has, and a 
variable set U G L initially empty. 

1: Cuts ^0 > Reset number of cuts 

2: for all / G L do 

3: 0(/) ^ > Reset number of other children 

4: end for 

5: for all c E S° do > Timewise 

6: li,l2^l{t) 

7: if [h G Uq] a [xb{li) = e] then 
8: Oih) ^ Oih) + 1 

9: if Oih) = O^ik) + 1 then 

10: Xdih) ^2 > Death time to new maximum 

11: Cuts ^ 1 

12: end if 

13: if Oih) > O^ih) then > h has too many other children 

14: li E U > We do not know the father of h 

15: end if 

16: else if [h E Uq] A [xbih) = then > Same, switching /i and h 

17: Oih) ^ Oih) + 1 

18: if Oih) = O^ih) + 1 then 

19: Xdih) ^ e 

20: Cuts 1 

21: end if 

22: if Oih) > O^ih) then 

23: h^U 
24: end if 

25: end if 
26: end for 

return Cuts, U, P 
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Conventionally, we may set V{B) to infinity: the border kills everything that 
hits it, and it will be easier to keep track in the algorithm. 

Lemma 3.2. With input satisfying the requirements given at the beginning 
of the section, the Algorithm 2 ends. 

Its output is a pretessellation satisfying the Late Events Property . More- 
over, each line has the same number of children as in the real tessellation 
T. 

Proof in appendix. 

The lemma states that the algorithm ends, but not that we have the real 
tessellation. I confess that I do not know whether we may have patholog- 
ical situations where there are several pretessellations with these numbers 
of other children and whose segments have the right endpoints to yield the 
same initialisation. However, we may circumvent the difficulty by carefully 
choosing the set Uq of lines whose parents we specify at input. 

For a given set f/o, we may write P or P{Uq) for the output pretessellation, 
S{Uo) the associated segments and s{Uo, I) the segment on line /. We also call 
D{Uq) G L the set of lines such that s{Uo,l) ^ s{l). We may also consider 
the sets Di,{Uo) and Dd{Uo) of lines with different birth and death times, 
that is (Z) 7^ x'[{l) {resp. x^{l) ^ Those two sets may intersect. 

Obviously D,{Uo) U D^iUo) = D{Uo) and D,{Uo) C Uo- 

Now, if we choose the right line and give its birth time, then there will 
be at least two less lines in Dh. that one and another. Namely we consider 
a line /i whose birth time was wrong, and look at its wrong parent I2. The 
line I3 whose birth time we correct is the one that kills this parent I2 in the 
real tessellation. Since the first line li had the wrong parent, its fake birth 
was later than this kill, so that the parent I2 now lacks a child. Hence there 
is another line U who will be born on I2 at the end of the algorithm. Since 
the death time of I2 is now right, it will really be the parent of U. Formally: 

Lemma 3.3. Let Uq, P and D{Uo), Db{Uo), Dd{Uo) defined as above. 
Then there is a line I such that 

#Z}fc(f/o\{/})<#/^6(f/o)-2 (14) 

Proof in appendix. 

An important remark is that if all the times of birth are right, then P 
is the real tessellation T. Indeed, since the times of birth are right and P 
satisfies the Late Events Property , the segments of P contain those of T. 
Since P is a pretessellation, segments do not cross, hence the times of death 
cannot be later than in the real tessellation T. 
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This remark will finally ensure that we may choose a labelling scheme 
with enough orphan lines to prove existence of our measures. We start by 
choosing as orphan lines the inner nodes of the tree of birth of odd or even 
generations, whichever the biggest. We then use the former lemma to obtain 
a working Uq. We get: 

Theorem 3.4. For any set of k lines, for any e > the number of T- 
tessellations built on them is at most: 



where C depends only on e. 
Proof in appendix. 

3.3 The T-tessellation model is well-defined 

Our combinatorial result implies existence of the Gibbs measure (??) in such 
cases. 

Theorem 3.5. Let H{T) be an energy on T such that H{T) > -Ci^L{T), 
for some real C and any tessellation T . Then for any intensity t, the Gibbs 
measure with density 



is well-defined and finite. 

Proof. We have to prove that the measure is finite. We denote by c any 
absolute constant. Using the bound (1), Stirling formula and Theorem 3.4, 
we get: 




(15) 



djw{T) = exp(-i7(T))dA{^(L(T)) 




oo 



i-Mk) exp (-r) jjAfik) 




< oo. 



□ 
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4 Optimality remarks and perspectives 



Though we have used very violent upper bounds at times, there is no way 
to get a substantially better combinatorial result. Indeed let us consider for 
some integer a < k the following k lines on a square domain [0, 1]^: 

y = ^— forAG[l,a] 
a + 1 

X = for A G [1, A; — al. 

k- a + 1 

How many different T-tessellations can we build on those lines? A lower 
bound is given by supposing that all horizontal segments are maximal, that 
is go from border to border. Then each of the vertical segments is between 
two consecutive horizontal lines, and hence of length l/(fc — a + 1). More 
significantly, this means each one can be at {k — a + 1) different places, 
independently from each other since the vertical lines do not cross. So that 
there are at least {k — a + ly different T-tessellations that can be built on 
those lines. If we take a = k — k/{\nk), we may conclude: 

Lemma 4.1. There are sets ofk lines such that the number of T-tessellations 
on those lines admits the following lower bound: 



U{k) > 



^ ^ fc-fc/(lnfc) 



In k 



If we want to get a better result and a tighter upper bound on the partition 
function, we then need to have a closer look on the usual topologies of the 
lines. That is an order of magnitude harder, but might be worth the effort. 
Indeed the previous worst-case example hinges heavily on having many lines 
crossing many segments, and topologically equivalent sets of lines have very 
low measure, looking like (/c2*'/c!)/(2A;)! of the space of all sets of k lines. 

By contrast, using very sloppy heuristics, we would expect that for most 
sets of k lines, the number of T-tessellations on those lines behaves like 

M=^rk. 

The idea is the following: let us take a segment away of the true tessellation. 
How many different segments may we put on the line to get a tessellation 
again? Neglecting problems of children and murders, this would be the num- 
ber of segments that the line cross, plus one. Now the probability of crossing 
a segment is essentially the length of this segment. So the number of crossed 
segments would be km, where m is the mean length of a segment. Now that 
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length is the interval between two successive segments a line cross, that is 
l/{km). So that m should be of order and for each new line, we have 

as many possibilities. 
Thus it seems likely that the method in this paper gives little information 
on the measure, except its very existence. 
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A Proofs for the algorithms 
A.l Proof of Lemma 3.1 

At the end of initialisation, we have the following properties: 

• Birth times are those of the tessellation for all lines /: (/) = xj{l). 

• Death times are overestimated for all lines /: x^{l) > 

• The number of murders M{1) for each line is that of the true tessellation. 

Indeed the first and third points are merely the input, and the death 
times are set to an upper bound at stage 1. 

What is important is that those properties will remain true throughout 
the loop that completes the algorithm, and this will yield by recurrence that 
at the end of the time e loop: 

• The remaining number of murders for each line M{1) is that of the true 
tessellation M'^{1, e) = # [m G L : a;J(m) = e(m, /) and e(m, /) > e}. 
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• Death times before e are right, that is: < e) ^ i^di^) ~ (0) 

We have to prove that if this is true before the c(/i,/2) loop, it will be 
true after it. 

Now, we enter the loop if and only if there is a death in the real tessella- 
tion. Indeed, in that case 

xi:{h) = xl{h)<c{h,h)<x^,{h)<x^{h), 

and the same for I2. If on the contrary there is no death, since segments 
do not cross (??), either xj{l) > e(/i,/2) for one of the two lines, and then 
this also holds for (/) = or one of the two lines / is already dead 

< e(/i,/2). Then x^{l) < by recurrence hypothesis. 

If we do not enter the loop, there are no changes to M or . On the 
other hand, there is no change to M-^, nor any new line whose death time is 
required to be right, so the conditions still hold. 

If we do enter the loop, then either = e(/i, I2), or a;J(/2) = ^2)- 

In the first case, using the recurrence hypothesis, M(/i) = M'^{li, e(/i, /2)) = 
0, and x^{li) is set to satisfying the second condition. The first 

condition is also still satisfied, since the only number of murders that changes 
for the real tessellation is that of I2, which decreases by one, since li is no 
more in the set of remaining murders. Symmetrically, if x^{l2) = c(/i,/2), 
then M(/i) = M'^ {li, z{li,l2)) > since it contains I2, and this number of 
remaining murders is decreased by one while x^{l2) is set to xJ(/2)- So that 
the recurrence hypothesis is transmitted. 

Since death times before e are right, after we hit the last event, all death 
times are right, that is x^{l) = for all lines /. Hence the output 

pretessellation is the real tessellation. 

A. 2 Proof of Lemma 3.2 

The proof is built on the fact that throughout the algorithm, after preinitial- 
isation, three conditions are fulfilled by every line / of the algorithm. The 
two conditions (7) and (8) of the Late Events Property , and the fact that if 
we are sure we do not know yet the parent of a line, its birth time is strictly 
later than its birth time in the real tessellation. In equation: 

{leU)^xai)>xUl). (16) 

We must then check "the three conditions" still hold each time we change 
a birth or death time or the set U. 
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First, let us check "the three conditions" hold after preinitialisation. 
Indeed on the one hand at stage (2) we set the death times of all lines to 
the maximum possible, that is the rightmost point of the domain. So that 

On the other hand, we know the birth times of the lines I not in Uq. For 
those (/) = xl{l). For the lines I ^ Uo, whose parent we do not know, we 
set their birth time to more than maximum possible at stage (5). So that 
(/) > xj{l). Since U is initialised as Uq, condition (16) is fulfilled. 

The "three conditions" still hold at the end of the whole initialisation. 
We shall even prove that the prototessellation at the end of initialisation is 
exactly the pretessellation described as P{Uq) in equations (10), (11) and 
(12). 

Indeed on the one hand the birth times are changed only at stages (10) 
and (12). 

Now these statements are not reached for the lines I ^ Uq, whose birth 
times we know: if t{l,m) = xf(m), then either / is really the parent of m 
and xj{m) = t{l,m) > xj^{l) = a;f(/), and we do not reach the change, or 
x^{m) has been changed at stage (10) or (12), at which time it was set to 
(/) = xj{l). Contradiction again. Hence times of birth of lines not in Uq 
are as described in equation (10) at the end of initialisation. 

For lines not in Uq, at stages (10) and (12), we set the birth time of a line 
/ G t/o to t{l,m) for some m ^ L. This happens once since afterwards the 
condition Xb{l) > t will never be satisfied. This does happen when m is a 
child of / in P. Since the first child of / is specified, according to requirement 
2, this is the event to which the birth time is set, and equation (11) is satisfied 
at the end of the initialisation. 

Since moreover m is the child of /, we know that x1{l) < t{l,m) = 
x^{l). In addition notice that the birth point is always included in the real 
tessellation segment s^(/): 

xr(/)<xi(o. (17) 

This condition also holds for the lines not in Uq since then xf (Z) = x^(/) < 
xj(/). It will hold throughout what's left of the algorithm, since birth times 
are never increased. In the latter case, the segment is then even included in 
the segment of the real tessellation T. 

On the other hand, death times are changed during initialisation only at 
stages (15), (17) and (23). We want to prove that at the end of initialisation, 
the death time is that mentioned when describing V , that is x^^^°'^ at equation 
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(12) . Since death and birth times can only be changed to e, the recurrence 
hypothesis imphes that the if condition Xb{li),Xb{l2) < £ < Xd{h),Xd{l2) is 
equivalent to it being true at the end of initialisation. Under that hypothesis, 
we enter the if statements with / G /(c) if and only if x^^^°\m) = e or 
x^^^°\l) = e. This happens exactly V{1) + 1 times, with the definition 

(13) . Each time save the last, V{1) is decreased by one, so that it really is 
at the last event, and x^{l) is set to x^^^°\l), transmitting the recurrence 
hypothesis and satisfying equation (12). 

Hence at the end of initialisation, the prototessellation P is really the 
pretessellation P{Uq). 

We now have to deal with the main loop, separated in functions Parent _seek 
and Cutting. ?? Not only do we have to check the "three conditions", but 
we shall also make use of the fact that P stays a pretessellation from now 
on: segments do not cross (3). 

Within Parent _seek, the changes occur at stages (4) and (5), and sym- 
metrically (8) and (9). Since any line whose birth time we change is excluded 
from U, condition (16) still holds. There are no change in death times, so 
that condition (8) still holds. 

If (7) still holds, P will stay a pretessellation. Indeed, the extended parts 
are then between times of birth and times of first child. These are included 
in the real tessellation T, so they cannot cross. A would-be crossing would 
then be with a segment s{m) that has not changed during the loop. But then 
the birth time would have been changed at e(/, m). Since a line has to be in 
U to have its birth time changed, and is excluded from U if it is changed, its 
birth time can change only once through one call of Parent_seek, finishing 
to show that segments still do not cross (3). 

We then have to ensure that birth times are still late (7). Since a birth 
time can change only once during one call and the events are tested reverse 
timewise, it is enough to prove that if / G ?7 at event xji^l), then (/) is 
set to xj{l). This is true, since the parent m of a line I in U G Uq is in 
L\Uq, thanks to requirement 2. So that this event is in the segment of the 
pretessellation: (m) = xj{m) < c{l,m) < xj(m) < x^{ni). We thus enter 
the if (3) or its symmetrical (7). 

The latter point also proves that when entering the function Cutting, all 
lines are born on a segment (4). 

Within Cutting, the changes happen at stages (17), (18), (23) and the 
symmetrical (17), (18), (23). No birth time is changed so condition (7) 
remains valid. To understand what is going on, let us consider an event 
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t{l,m) where I G Uq, m ^ L, and = c{l,m) at input. Then since the 

birth time is included in the real segment (4), either / is really a child of m, 
that is xj{l) = e(/, m), or since moreover segments in T do not cross (3), the 
segment m dies before: xj(m) < e(Z, m). Now we change death time at stage 
(10) (or (17)). The counter O is the number of [/-children of m at input 
that are born at e(/, m) at the latest. So that when we change its death time 
(10), there are 0^(m) such children strictly before e(/,m). Since it cannot 
have more children, we obtain xj(m) < e(/,m) and the condition on death 
times (8) is still fulfilled. Besides, when we add a line / to U, at stages (14) 
and (23), we still have (/) < x^(/), fulfilling condition (16): indeed this 
happens at points e(/,m) where 0{m) > 0^{m), so that xj(m) < t{l,m) 
thanks to the former remark, and m cannot be the parent of I. 

Notice that 0{B) > 0^{B) is impossible thanks to the late births (7). 

Since segments are only shortened during Cutting, they will not cross 
and P stays a pretessellation. 

Finally notice that Cuts is set to one, at stages (11) or (18), if and only 
if at least one line's death time gets earlier. 

This remark implies that the algorithm ends: indeed death times can only 
decrease, and they take values in a finite set, so there will be an iteration 
when they do not change, there are no cuts. The algorithm ends at this 
point. 

What can we say about the final pretessellation returned by the algo- 
rithm? 

To begin with, it is really a pretessellation that fulfills the Late Events 
Property . We have followed the conditions throughout the algorithm. 

Moreover, each line has the same number of children as in the true tessel- 
lation. Indeed the algorithm ends only when no cuts are made in the Cutting 
function, that is when all lines have at most as many children as in the true 
tessellation. The other side of the inequality comes from the fact that the 
sum of other children in that function is bounded from below by the sum of 
other children in the true tessellation: each line is born on another line or on 
the border, and since birth times are overestimated, they cannot mistakenly 
be thought to be born on the border. 

A. 3 Proof of Lemma 3.3 

A first remark is that for any subset Ui of Uq, the Late Events Property holds 
for P{Uo) with respect to P{Ui). Indeed there would be no change when 
running Algorithm 2 with Uq as input if P{Ui) was the real tessellation. So 
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that x^^^°\l) > x^^^^\l) > xj{l) for all lines. This implies that Dh{Ui) is a 
subset of Dh{Uo). 

Another remark, already made in the proof of the algorithm, is that when 
= e(Z,m), either m is the real parent of /, or a;J(m) < e(/,m). This 
implies that if xj(m) = x^{m), all its children in P are real children. Since 
moreover it has the same number of children as in T, its children are exactly 
the right children. 

We now consider li G DbiJJo) and its fake parent Z2, so that (/i) = l2- 
Then t^li.h) > x^{l2)- Moreover t^hjs) = 3^J(^2) and I3 G D^iUo). Indeed 
Xb{h) > ^{h,h)- We take I3 as the line I in the lemma. Since ^3 ^ ^^\{^3}, 
we know that I3 ^ Dh{Uo\ {h})- 

Now ^3 is no longer a child of I2. But the number of children of I2 at the 
end of the algorithm is fixed, equal to that in the real tessellation T. So that 
there is a line U that is a child of I2 in P{{Uo\ {h}) and was not in P{Uq). 
Since moreover I2 is killed by the right line in P{Uq\{1^}), we know that 
/4e A(f/o)and h^D,{Uo\{h}). 

Finally, since Di,{Uq\ {h}) £ Dh{Uo), we may write: 

#D,(f/o\ {/}) < i^D,{Uo) - # {/3, ^4} = #A(t/o) - 2. 
A. 4 Proof of Theorem 3.4 

We may rebuild any tessellation T on the lines L with Algorithm 2 and the 
input given in the "requirements" mentioned at the start section 2. However 
we have to make it more specific, namely to describe Uq the set of orphan 
lines. 

This set is obtained by looking at the tree of births. On the one hand, 
we must specify the parent of each leaf, since we cannot give their first 
child. We add to this the interior nodes either of even, or of odd generations, 
whichever the smaller. Thus the lines f/j without specified parents are at 
most Ui = {k — z)/2, where z is the number of leaves. This is a valid labelling 
as the lines whose parents are not specified have all their children marked as 
such, since the parity of their generation is different. Hence requirement 2 
for the input is satisfied. 

Now we run the algorithm. We may not find the true tessellation, and 
have a set of lines Dh{Ui) whose parent is wrong. This set is included in 
Ui. We then remove a line from Uq as in Lemma 3.3 and run the algorithm 
again. And we iterate until we obtain the true tessellation. Since Dh{Uo) is 
at least two elements smaller at each steps, we have to remove at most u/2 
lines from our initial Ui to get a set Uq we may use as input in Algorithm 2 
to obtain the true tessellation. So that its cardinal is at least u = {k — z)/4:. 
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Let's now go over each element of the input in order. 

First labelhng each hne as "leaf" or "not leaf". As already mentioned, this 
yields at most 2^ possibilities. 

Second labelling each line not in t/o with its parent. Since a leaf cannot 
be a parent, there are at most {k^\ — z)^~'^ < (k + 1 — z)^^^'^^'^^^ possibilities. 

Third labelling each line with its number of virtual murders yield at most 
4*^ possibilities, as already mentioned. 

Fourth labelling each line with its number of other children yield at most 

possibilities, as already mentioned. 

Thus we may give the following upper bound on the number of different 
T-tessellations on k given lines, using C for any absolute constant: 



where we have used the following bound on the supremum in the right-hand 
side: take the derivative in z of the logarithm, and we see that the maximum 
is attained when 



For big k, this implies k/\n{ky "^yi + k — z>4k/ ln{k). We then replace 
by the right bounds in the exponent and the basis. 



ATik) < C'- sup {k + 1- 




l<z<k 
/ U \ fc-fc/(lnfc) 



{l + k-z){l + ln(l + k-z))=Ak-l. 
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